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We present a calculation of current-induced torques in metallic magnetic multilayers derived from 
the spin-dependent transmission and reflection properties of the magnetic layers. A scattering 
formalism is employed to calculate the torques in a magnetic-nonmagnetic-magnetic trilayer, for 
currents perpendicular to the layers, in both the ballistic and diffusive regimes. 
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I. INTRODUCTION 

Stacks of alternating ferromagnetic and nonmagnetic 
metal layers exhibit giant magnetoresistance (GMR), 
because their electrical resistance depends strongly on 
whether the moments of adjacent magnetic layers are 
parallel or antiparallel. This effect has allowed the devel- 
opment of new kinds of field-sensing and magnetic mem- 
ory devices. tl The cause of the GMR effect is that con- 
duction electrons are scattered more strongly by a mag- 
netic layer when their spins lie antiparallel to the layer's 
magnetic moment than when their spins are parallel to 
the moment. Devices with moments in adjacent mag- 
netic layers aligned antiparallel thus have a larger overall 
resistance than when the moments are aligned parallel, 
giving rise to GMR. This paper discusses the converse 
effect: just as the orientations of magnetic moments can 
affect the flow of electrons, then by Newton's third law, 
a polarized electron current scattering from a magnetic 
layer can have a reciprocal effect on the moment of the 
layer. As proposed by BergerB and Slonczewskijj an elec- 
tric current passing perpendicularly through a magnetic 
multilayer may exert a torque on the moments of the 
magnetic Jayers. This effect which is known as "spin 
transfer" ,□ may, at sufficiently high current densities, al- 
ter the magnetization state. It is a separate mechanism 
from the effects of current induced magnetic fields. Ex- 
perimentally, spin-ouigent-induced magnetic excitations, 
such as spin- waves pa and stable magnetic reversal,!™ 
have been observed in multilayers, for current densities 
greater than 10 7 A/cm 2 . 

The spin-transfer effect offers the promise of new kinds 
of magnetic devices,B and serves as a new means to ex- 
cite and to probe the [dynamics of magnetic moments 
at the nanometer scale.E3 In order to controllably utilize 
these effects, however, it is necessary to achieve a bet- 
ter quantitative understanding of these current-induced 
torques. Slonczewski has presented a derivation of spin- 
transfer torques using a l—D WKB approximation with 
spin-dependent potentials,ubut his calculations only take 
into account electrons which are either completely trans- 
mitted or completely reflected by the magnetic layers. 
For real materials the degree to which an electron is 
transmitted through a magnetic/nonmagnetic interface 



depends sensitively on the-|p,atching of the band struc- 
tures across the interface. tilE3 It is the goal of this paper 
to incorporate such band structure information together 
with the effect of multiple reflections between the ferro- 
magnetic layers, into a more quantitative theory of the 
torques generated by spin-transfer. TMs could be done 
using the formalism of Brataas et oZ,c3 which is based 
on kinetic equations for spin currents. Instead we choose 
to employ a modified Landauer-Buttiker formalism, in 
which we model the ferromagnetic layers as generalized 
spin-dependent scatterers. The calculations are carried 
out for a quasi-one dimensional geometry, for which we 
derive formulas for the torque generated on the magnetic 
layers when a current is applied to the system, for cither 
ballistic or diffusive non-magnetic layers. The main dif- 
ference between our approach and Ref. [l3] is that in our 
case, scattering in the normal layer is phase coherent, 
whereas Ref. O assumes phase relaxation. However, in 
the case of a diffusive normal metal layer and for a large 
number of transverse modes, the two approaches would 
give the same answer. 

The paper is organized as follows: in Section ||, we 
present an intuitive picture (adapted frorrrra) of how 
spin-dependent scattering of a spin-polarized current pro- 



duces a torque on a magnetic element. Section III is de 



voted to the introduction of the scattering matrix formal- 
ism for the spin-flux. This formalism is then used in Sec- 
tion IV to calculate the torque in a Ferromagnet-Normal- 



Ferromagnet (FNF) system where the normal part is dis- 
ordered (diffusive). Section ^ contain a discussion of the 
results. Details of our calculation are presented in Ap- 
pendix A. In appendix B, we derive the torque for an 
FNF system where transport in the normal layer is bal- 
listic, rather than diffusive. 



II. PHYSICAL IDEA 

In this section we will present a simple intuitive picture 
of the physics behind the spin-transfer effect. The con- 
nection between current-induced spin-transfer torques 
and the spin-dependent scattering that occurs when elec- 
trons pass through a magnetic/nonmagnetic interface can 
be illustrated most simply by considering the case of a 
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Torque ~ sin 6 



FIG. 1. Schematic of exchange torque generated by 
spin-filtering. Spin-polarized electrons are incident perpen- 
dicularly on a thin ideal ferromagnetic layer. Spin-filtering 
removes the component of spin angular momentum perpen- 
dicular to the layer moments from the current; this is ab- 
sorbed by the moments themselves, generating an effective 
torque on the layer moments. 



spin-polarized current incident perpendicularly on a sin- 
gle thin ferromagnetic layer F, as shown in Figure 1. The 
layer lies in the y — z plane, with its magnetic moment 
uniformly pointed in the +z direction, and we assume 
that the current is spin-polarized in the z — x plane at an 
angle 9 to the layer moments. The incoming electrons can 
therefore be considered as a coherent linear superposition 
of basis states with spin in the +z direction (amplitude 
cos(0/2)) and -z direction (amplitude sin(#/2)). For this 
initial discussion we will assume that the layer is a perfect 
spin-filter, so that spins aligned with the layer moments 
are completely transmitted through the layer, while spins 
aligned antiparallel to the layer moment are completely 
reflected. For incident spins polarized at an angle 9, the 
average outgoing current will have the relative weights 
cos 2 (9/2) polarized in the +z direction and transmitted 
to the right and sin 2 ((9/2) polarized in the — z direction 
and reflected to the left. Consequently, both of the out- 
going electron spin fluxes (transmitted and reflected) lie 
along the z axis, while the incoming (incident) electron 
flux has a component perpendicular to the magnetiza- 
tion, along the x axis, with magnitude proportional to 
sin 9. This ^-component of angular momentum must be 
absorbed by the layer in the process of filtering the spins. 
Because the spin-filtering is ultimately governed by the 
s — d exchange interaction between the conduction elec- 
trons and the magnetic moments of the layer, the angular 
momentum is imparted to the layer moments aad pro- 
duces a torque on them. This exchange torque,li-3 which 
is proportional to the electron current through the layer 
and to sin (9, is in the direction to align the moments with 
the polarization of the incident spin current. 

The symmetry of this model precludes any generation 
of torque from the spin-filtering of a current of unpo- 
larized electrons. To generate the effect, then, a second 
ferromagnetic layer is needed to first spin-polarize the 



current, see Fig. |2[ In that case, spin angular momen- 
tum is transferred from one layer to the flowing electrons 
and then from the electrons to the second layer. However, 
the torques on the two layers are not equal and opposite, 
as spin angular momentum carried by the electrons can 
also flow away from the layers to infinity, see Fig. 0. 
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FIG. 2. Qualitative picture of asymmetry of spin-transfer 
torque with respect to current bias in a FNF junction. For 
left-going electrons (2a), initially polarized by a magnetic 
layer F a , the moments of layer Ft experience a torque so as 
to align them with layer F a . The electron current reflected 
from layer Ff,, in turn, exerts a torque on layer F a so as to an- 
tialign it with the moment of layer Fj,. Subsequent reflections 
between the layers reduce but do not eliminate this torque. 
If the current is reversed (26) the overall sign of the torque 
is reversed, encouraging the moment of layer Fj, to align an- 
tiparallel with layer F a . 

The presence of this second layer has the additional 
effect of allowing for multiple scattering of the electrons 
between the two layers, which gives rise to an explicit 
asymmetry of the torque with respect to current direc- 
tion. This asymmetry is an important signature which 
can be used to distinguish spin-transfer-induced torques 
from the torques produced by current-generated mag- 
netic fields. To see how the asymmetry arises, consider 
the ferromagnet-normal-metal-ferromagnet (FNF) junc- 
tion shown in Figure |2| It consists of two ferromagnetic 
layers, F a and Fb, with moments pointing in directions 
rh a and mj,, separated by a normal metal spacer N. Nor- 
mal metal leads on either side of the trilayer inject an 
initially unpolarized current into the system. When the 
current enters the sample from the left (Fig. |]a), elec- 
trons transmitted through F a will be polarized along rh a . 
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As long as the normal metal spacer is smaller than the 
spin-diffusion length (100 nm for Cu), this current will 
remain spin-polarized when it impinges on Ft, and will 
exert a torque on the moment of F& in a direction so as 
to align rhb with rh a . Repeating the argument for Ft,, wc 
find that the spin of the electrons reflected from layer Ft, 
is aligned antiparallel to rhb, and, hence in turn, exerts a 
torque on the moment of F a trying to align rh a antiparal- 
lel with rhb- (Subsequent multiple reflections of electrons 
between F Q and Ff, can serve to reduce the magnitudes of 
the initial torques, but they do not eliminate or reverse 
them, as the electron flux is reduced upon each reflec- 
tion.) When the current is injected from the right, the 
directions of the torques are reversed: Now the flow of 
electrons exerts a torque on F a trying to align its mo- 
ment parallel with rhb, while it exerts a torque on Ft, so 
as to force the moment in layer Ff, antiparallel with rh a . 

In the remainder of the paper, we assume that rhb 
points in the +z direction, while rh a differs by a small 
angle in the x — z plane. (For thin films, demagnetizing 
forces will in general cause the y — z plane to be preferred, 
but this produces no change in our argument. We present 
the case that is easier to draw.) The overall effect of a 
left-going flow of electrons then, is to exert a torque ft, on 
Ft, in the +x direction. If we reverse the current, so that 
electrons pass through Ff, first (Fig. 2b), the torque on 
Ff, is only exerted by the electrons after they have been 
reflected from F a . As seen before, the electrons reflected 
from F a have polarizations opposite to rh ai so that the 
torque on Ff, is in the —x direction. 

In Refs. |^ and |?], the layer F a was taken to be 
much thicker than Ft,, so that intralayer exchange and 
anisotropy forces will hold the orientation of rh a fixed. 
In that case, one is only interested in the torque on Ft,, 
which serves to align rhb either parallel or antiparallel 
with the fixed moment rh a depending on the current di- 
rection. This asymmetric current response has been em- 
ployed in both a point-contact geometryu and in a thin- 
film pillar geometryQ to switch the moments in FNF tri- 
layers from a parallel to an antiparallel configuration by a 
current pulse in one direction, and then from antiparallel 
to parallel by a reversed current. For weakly-interacting 
layers, either orientation can be stable in the absence of 
an applied current, so that the resistance versus current 
characteristic is hysteretic, and the devices can function 
as simple current-controlled memory elements. 



III. SPIN FLUX AND TORQUE IN THE 
SCATTERING APPROACH 

Treating the ferromagnetic layers as perfect spin filters 
provides important qualitative insights into spin-transfer, 
but for a complete qualitative and quantitative picture, a 
more general approach is required. In this section, we in- 
troduce a scattering matrix description of the FNF junc- 
tion which allows us to deal with non-ideal (magnetic and 



non-magnetic) layers. Our goal is to relate the torque ff, 
exerted on layer Ff, by an unpolarized incident electron 
beam to the scattering properties of the layers. Although 
we shall restrict our formulas to the FNF junction (see 
Fig. ^]), our method is applicable for an arbitrary array 
of magnetic-nonmagnetic layers. 

We first introduce the spin flux J in the a;-direction 
(the direction of current flow): 



J(x) = - — Im 
v ' 2m 



dydz 



0t( 



,]a Tx^ x) 



(i) 



where 4>{x) is a spinor wavefunction and a the vector of 
Pauli matrices, 



<f>(x) 




Note that although Eq. (|l|) bears close formal analogy to 
the particle current, no local equation of conservation can 
be written for the spin flux, since in general, the Hamil- 
tonian docs not conserve spin. Specifically, the magnetic 
layers can act as sources and sinks of spin flux, so that the 
spin flux on different sides of a F layer can be different. 




FIG. 3. Schematic of the setup used for the definition of 
the scattering matrices of the F and N layers. The two layers 
F a and Ft are ferromagnetic layers whose magnetic moment 
is oriented as shown in the bottom of the figure. The layer N 
is a nonmagnetic metal spacer. Amplitude of left and right 
moving propagating waves are defined in fictitious ideal leads 
0, 1, 2 and 3 between the layers and between the layers and 
the reservoirs. 



A. Definition of the scattering matrices 

Fig. U shows the FNF junction where (fictitious) per- 
fect leads (labeled 0, 1, 2 and 3) have been added in be- 
tween the layers F and N and between the F layers and 
the electron reservoirs on either side of the sample. The 
introduction of these leads allows for a description of the 
system using scattering matrices. In the perfect leads, 
the transverse degree of freedom are quantized, giving N 
propagating modes at the Fermi level, where N ~ A/X F 
A being the cross section area of the junction and Xf 
the Fermi wave length. Expanding the electronic wave 
function in these modes, we can describe the system in 
terms of the projection ^^l/r of the wave function onto 
the left (right) going modes in region i. The ^i^L/R are 
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2iV-component vectors, counting the TV transverse modes 
and spin. The amplitudes of the wave function in two 
neighboring ideal leads are connected through the scat- 
tering matrices Sb, S a and Sn, that relate amplitudes 
of outgoing modes and incoming modes at the layer (see 
for example Ref. [lB] for a review of the scattering matrix 
approach) , 



*3L 
®2R 

*2L 



Sb 

S a 



*3fi 
®IR 
^2R 



(2a) 
(2b) 
(2c) 



The scattering matrices 5(, , S a and S*n are 47V x 4TV uni- 
tary matrices. We decompose Sb into 27V x 27V reflection 
and transmission matrices, 



S h = 



n t' b 
tb ri 



(3) 



with similar decompositions of S a and Sn- Normalization 
is done in such a way that each mode carries unit current. 
Due to the spin degree of freedom, the reflection and 
transmission matrices can be written in terms of four 
TV x TV blocks: 



r = I r bn r bU 

nu nn 



(4) 



where the subscripts |, J. refer to spin up and down in the 
z-axis basis. 

The scattering matrix of the magnetic layers depends 
on the angle the moments may make with the z-axis. 
The matrix S a (0) is related to S a (9 = 0) through a rota- 
tion in spin space: 

r«(0) = Re r a (0) R-e, r' a (6) = R e r' a (0) R^g, 
ta(0) = Re t o (0) R-e, t' a (9) = R e t' a (0) R-e (5) 



where 



Re — 



cos ■ 
sin ' 



1 



N- 



(6) 



The non-magnetic metallic layer will not affect the spin 

states, i.e., r-Nfj. = ?"nj,t = an d r NTT = r Nj.i- 

We need to keep track of the amplitudes within the 
system in order to calculate the net spin flux deposited 
into each magnetic layer. Therefore, we define 2TV x 2TV 
matrices Tf' R and Af' R (i = 0, 1,2,3) so that we may 
express all the ^^l/r as a function of the amplitudes 
incident from the two electrodes (regions and 3): 



T iR A iR J { $> 3R 



(7) 



with the convention that Tql — A3/? — 1 and T3H = 
Aql = 0. In order to calculate the torque exercised on 



layer for a current entering from the left, we need the 
matrix, T2L- To simplify the notations in the rest of the 
paper, we write 



ci = r 



2L- 



(8) 



The matrix fl relates VE^L to the incoming amplitudes 
^ol coming from the right. To calculate it, we put = 
0, then, using Eq.(|^), we get the equations: 



^XR = t n ^ 2 R 
= r n ^ 2 R 

%2R 
*3L 



t' b *2L, 



from which we obtain: 
1 



1 



-t' 



'n'b 



1 - r n t n r 



■ , 1 H — r r' 
a°ni b l—r n r'. « a n 



(9) 



(10) 



B. Spin flux response 

Let us now connect our system to two unpolarized elec- 
tron reservoirs on its two sides, as shown in Fig. In 
equilibrium, the modes in the reservoirs are filled up to 
the fermi level tF- We want to calculate the spin current 
that is generated when the chemical potential in the left 
(right) reservoir is slightly increased by 6/13 (<fyto)- The 
spin current Ji is the difference of the left going and right 
going contributions. For each of the region i = 0,1,2,3, 
we find from Eq.(EI) and Eq.(0): 



d-k = J_ Re 
dfi3 An 



TroT^r^ - TroTi/,r 



and 



dJj _ 1 
djio An 



Re 



TraA iR Al R - TraA^A^ 



(11) 



(12) 



Derivation of Eq. ( |ll[ ) and Eq. (|12|) proceeds analogously 
to the derivation of the Landaucr formula for the 
conductance.cJ 
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FIG. 4. The FNF junction is connected to two reservoirs. 
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C. Torque exercised on layer Ft, 

If the spin flux on both sides of F& (region 2 and 3) is 
different, then angular momentum has been deposited in 
the layer F&. This creates a torque t% on the moment of 
the ferromagnet, 



n = J3- Ji- 

Setting 5(j,q = — eVo, we have: 

dfb __ 
dV ~ 4tt 



-Re Tr 



2N 



with 



t'Jat' h 



't - ' 



(13) 



(14) 



(15) 



This equation can be simplified further if the spin- 
transfer effect is. due entirely to spin-filtering (as argued 
by SlonczewskitJ) as opposed to spin-flip scattering of 
electrons from the magnetic layers. That is, if we as- 
sume that r bn = r bU = r aU (6 = 0) = r aiT (6» = 0) = 0, 
then: 



dV Q 



2tt 



-Re Tr 



JV 



(n TT n+ +n n n] [,) 



(1 



r &TT r bll 



tl 



f T ) 



(16) 



We will comment briefly in the conclusion of this pa- 
per about the physical implications of including the off- 
diagonal spin-flip terms, as well. 

We note that, as there is no spin flux conservation 
in this system, dJi/dfi^ can be different from —dJi/dfio 
and, hence, there can be a non zero spin flux even when 
the chemical potentials are identical in the two reser- 
voirs. The existence of a zero-bias spin-flux and the re- 
sulting torques reflect the well-known itinerant-electron- 
mediated exchange interaction (a.k.a. the RKKY inter- 
action) between two ferromagnetic films separated by a 
normal-metal spacer. This interaction cap,ipJact be un- 
derstood within a scattering framework. EZTlJ The zero- 
bias torque has to be added to the finite-bias contribution 
(given by Eq.(|l6|)). Since the former is typically a factor 
TV -1 smaller and vanishes upon ensemble averaging (see 
section IV and Ref. |l5|), we henceforth neglect the zero- 
bias contribution to the torque and restrict our attention 
to the bias induced torque, for which we have 

df b _ df b 
dV Q ~ 8V 3 

up to a correction of order N^ 1 . 



IV. AVERAGING OVER THE NORMAL LAYER 



Via Eq.(p_6[), the torque on the moments of the fer- 
romagnetic layers F a and F& not only depends on the 



scattering matrices S a and Sb of these layers, but also 
of the scattering matrix Sn of the normal metal layer 
in between. If the normal layer is disordered, f a and f b 
depend on the location of the impurities; if N is ballis- 
tic the torque depend sensitively on the electronic phase 
shift accumulated in N. In general, sample to sample fluc- 
tuations of the torque will be a factor N^ 1 smaller than 
the average.E3 Hence, if N is large (N > 10 3 in the ex- 
periments of Ref. [?]) , the torque is well characterized by 
its average. In this section, we average over Sj$ for the 
case where N is disordered. The case of ballistic N is 
addressed in appendix B. After averaging, the zero-bias 
spin transfer current, corresponding to the RKKY in- 
teraction described above, vanishes, and only the torque 
caused by the electron current remains. Because all ef- 
fects of quantum interference in the N layer will disappear 
in the process of averaging (to leading order in the num- 
ber of modes N), the results we derive are unchanged 
if the reflection and transmission matrices include pro- 
cesses in which the energy of the electron changes during 
scattering JL 3 ] in addition to the elastic proceses normally 
considered in scattering matrix calculations. 



A. Averaged torque 

The scattering matrix of the normal layer pan be writ- 
ten using the standard polar decomposition:^] 



Sn 



u 
V 



U' 
V 



(17) 



where U, V, U' and V are 2N x 2N unitary matrices and 
T is a diagonal matrix containing the eigenvalues of tnt^. 
Since is diagonal in spin space, we find that U, U' , V 
and V are block diagonal. 



U 



u 
u 



U' 



u' 
u' 



(18) 



and similar definitions for v and v' . In the isotropic 
approxi m at i o nl 15 H 2 H , the N x N unitary matrices u,u',v 
and v' are uniformly distributed in the group U(N). (The 
outer matrices in 

Eq.(0) thus mix the modes in a ergodic 
way while the central matrix contains the transmission 
properties of the layer, which determine the average con- 
ductance of N.) 

We want to average Eq. (|l^) over both the unitary ma- 
trices and T. A diagrammatic technique for such aver- 
ages has already been developed in Ref. ^2] and can be 
used to calculate (dfb/dVo) in leading order in 1/N. It is 
a general property of such averages that the fluctuations 
are a factor of order N smaller than the average. This 
justifies our statement above, that the ensemble averaged 
torque is sufficient to characterize the torque exerted on 
a single sample. Details of the calculation are presented 
in Appendix A. 

The resulting expression for (dfb/dVo) can be written 
in a form very similar to the one for Eq. (|l6|) if one uses 
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a notation that involves 4x4 matrices. To be specific, 
to each 2N x 2N matrix A appearing in Eq.(16) and 
Eq. (|Io|) , we assign a 4 x 4 matrix A as, 



A = —Tr N [A ® A*] , 



(19) 



where Tr^v means that the trace has been taken in each 
the N x N blocks, or in extenso: 



A=-Tr N 



/ A TT Aj T A n A^ A Ti Aj T A n A^ \ 
A TT A ? A tt aJ x A n A T A U A* 
A^AL A iT A T n A U AL A n Aj l 

v%4 T %4i ^u4t a u4iJ 



We also define E by, 



(20) 



t' 



t' b n. 



(27) 



We would like to note that, while our theory started 
from a fully phase coherent description of the FNF tri- 
layer, including the full 4N x AN scattering matrices of 
the FN interfaces, the final result can be formulated in 
term of 2 x 4 parameters, represented by the matrices f a 
and f' b ( 2 x 16 parameters in case of spin- flip scatter- 
ing) . Such a reduction of the number of,degrees of free- 
dom was also found by Brataas et a/.,EJ although their 
starting point is an hybrid ferromagnetic-normal metal 
circuit with incoherent nodes. This confirms the state- 
ment at the beginning of this section, that for a diffusive 
normal-metal SDacer all effects of quantum interferences 
are washed out.ta The difference between our approach 
and the one of Ref. [r] is important in the case of the 
ballistic normal layer, see Appendix B. 



E = Tr 



A' 



/ E^f 








2n 



(21) 



The average over the transmission eigenvalues T follows 
if we note that the average of a function is the function of 
the average, to leading order in l/iVE3 Thus the average 
over T amounts to the replacement 



5n 
N 



I4 and f„ 



1 



fjN 

N 



(22) 



where <?n is the conductance of the normal layer and I4 
is the 4x4 unit matrix. Using these "hat" matrices, the 
result has now the simple form: 



Wo' 



47T 



-Re Tr 4 



EO 



where (compare to Eq.(10)), 
1 



n = - — — t' n — . - j - — . , 

b 1 — r a t n r b i_f n fi t n — r a r n 



1 - r n r 



-t 



(23) 



(24) 



Equation ( p3| ) is the main result of this paper. In the 
absence of spin-flip scattering, it reduces to 



(S)=-^Re (Si:,, <»,,. 

OVq Z7T 



x Trjv(l -r' m r 



't 



E 6TT Hi' 



(25) 



The same formalism can be used to calculate the conduc- 
tance g of the system using the Landauer formula. One 
gets: 



(<?> = 



Ne 2 



(26) 



B. Symmetries 

Before we proceed with a further analysis of Eq. (^3|) , 
we identify the different symmetries of the torque. Due to 
the conservation of current, the total torque deposited on 
the full system is anti-symmetric with respect to current 
direction: 



dr b dTo_ 
dV dV 



On 
dV 3 



dV 3 



(28) 



Eq.(p8|) holds before averaging. But, as pointed out in 
section [II, equality for each of the torques f a and t% 



separately holds only after averaging, 



Wo' 



y dv 3 > 



(29) 



sample to sample fluctuations of drb/dVo and —dfb/dV 3 
of relative order 1/N are in general different. Thus, for 
N 3> 1, our calculation can be used to compute the linear 
response of the torque to a small bias voltage: 



n 



1 dn 

Wo 



>(Vo W 3 ). 



(30) 



In our geometry, where F a and F;, are ini-z plane, the 
only non-zero component of the torque is r b . The torque 
vanishes when the moments are completely aligned or 
antialigned (all the matrices are diagonal in spin space 
and therefore no x-component of the spin can be found) . 
Around these two limits, the torque is symmetric in re- 
spect to the angle (6-^—6 and it — 9 — > tt + 9). There is 
no symmetry between 9 and tt — 9. In addition, the two 
layers are not equivalent and exchanging the scattering 
matrices of F a and F b also changes the torque. 



t' being the total transmission matrix: 



G 



C. Discussion of some limiting cases 

Eq.(p3|) can be simplified in some particular cases. Let 
us start with the simplest case of ideal spin filters, so 
that majority (minority) spins are totally transmitted 
(reflected) by either layer. Equation (^3|) then reduces 
to 



W 



e c/n sin 9 
Att 3 + cos 9 



h , . tan 0/2 , „ 
"i^^T"' (31) 



where (g) is the average magnetoconductance, cf Eq.(26) 
- 2 4 cos 2 (9/2 



(9) 



T ffN "3 



cos f 



(32) 



Equation (|3l| ) reproduces a result of Slonczewski.H As 
expected, for left-going electrons (Vq < 0) the torque 
is positive, so it acts to align the moments of the two 
magnetic layers, see section O. 

Let us now consider the case of weak s — d exchange 
coupling, i.e., when the scattering coefficients depend 
only weakly on spin. We continue to assume that no spin- 
flip scattering occurs in the ferromagnetic layers. We de- 
fine g a and gb as the average conductance per spin of the 
two layers (in unit of e 2 /h). Then, the conductance of F a 
alone is g a + 5g a and g a — Sg a for respectively the major- 
ity and minority spins, which defines the spin scattering 
asymmetry Sg a . In that case, we get to lowest nontrivial 
order in 5g a and Sgb'- 



2e 2 

= IT v 



g^g a gb 



•5N(5a + 5b-2^f) 



g^SgaSgb cos ( 



( 2 #+3N(ffa 



(33) 



and exit the system from the right hand side. There- 
fore these electrons cross the normal layer at least twice 
and this leads to the factor g^. On the other hand the 
conductance is linear in ^n- Therefore in order to maxi- 
mize the torque deposited per current, one has to use the 
cleanest possible normal metal spacer. (This statement 
is true in this limit of weak filtering, but not in general, 
see section |y|.) Note that in the previous case (perfect 
spin-filtering) the torque is proportional to <7n instead of 
the expected g^. Indeed, in that case, once the electron 
has been reflected by the layer Fb, it cannot go through 
Fb which works as a perfect wall for it. Therefore cur- 
rent conservation implies that it goes out of the system 
through the right. For <?n <C N, the torque is actually 
proportional to g 2 ^ for arbitrary spin asymmetry (except 
perfect filtering), and one gets: 



chf 
Wo 



) oc g^sin9, g N < N, 



(35) 



the factor of proportionality being a complicated function 
of the transmission probabilities of the layers. 



V. APPLICATION TO CURRENT-DRIVEN 
SWITCHING OF MAGNETIC DOMAINS 

In this section, we consider the general solution Eq. (|22 
for the spin-transfer torque. We first address strongly 
polarizing systems and then calculate torques for scatter- 
ing parameters more appropriate for the transition metal 
trilayers that can be studied experimentally. As the tri- 
layer devices are primarily current-driven, we calculate 
the torque per unit of current /, 

7? ,19t x , 



and 



x gdv ' 



Ab\ 

Wo' 



g^g a 5g 2 smt 



2tt2(1 



5N(ffa + .96 



JV 



(34) 



This last formula shows that: 

(i) The torque is not symmetric with respect to inter- 
changing the layer F a and Fb, in contrast to the conduc- 
tance. If one changes Sg a to —5g a , the sign of the torque 
is reversed. However, (dr^/dVo} oc Sg 2 , so if one changes 
5gb to —Sgb, the sign of the torque is unchanged. The 
sign of the torque on a ferromagnetic layer therefore de- 
pends on whether the other layer is a positive or negative 
polarizer, but not on the sign of filtering for the layer ex- 
periencing the torque. We have verified that this is true 
also in the general case. This point explains why the two 
layers can not be treated on an equal footing. 

(ii) We see that <?n appears through its square. In- 
deed, in order for some spin to be deposited in the layer 
Fb, some left going electrons have to be reflected by Fb 



The torque is measured in units of 

Eq. ( pl|) of the previous section gives the torque per 
unit current for the case that both layers F a and F& are 
perfect polarizers. The main feature of this system is 
that the 9 dependence of the torque is not of a simple 
sin 9 form, and that the torque per unit current diverges 
at 9 = 0. In Fig. ||, we look at what happens when one 
of the layers (Fb) is a nearly perfect polarizer while the 
other one is not. Although the divergence at 9 = tt is 
regularized, t^/I remains sharply peaked near 9 = tt. 
This is relevant for the critical current needed to switch 
the magnetization of Fb from 9 = tt to 9 = 0. Recall that 
the switching of the domains follows from a competition 
between the spin-transfer torque on the one hand and 
restoring forces from local fields, anisotropy, exchange 
coupling etc. (The competition between these forces has 
been considered phenomenologically in Ref. |§|,p3| using 
a phenomenological Landau-Lifschitz-Gilbert equation.) 
The torques for 9 close to and tt determine the critical 
currents to overturn a metastable parallel (antiparallel) 
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FIG. 5. Torque per unit current for the case where Fj, is 
a nearly perfect polarizer (|if,-]-| 2 = 0.999, |tb|| 2 = 0.001) and 
F a is not (|t aT | 2 = 0.3 |i a || 2 = 0.01). (solid line). The dashed 
line shows the case of perfect polarizers, see Eq.(|3~ij). Torque 
in measured in unit of 

27re 



alignment of the moment in F a and F;,. Hence the critical 
current should be different at 9 — and 9 — ir. 

In Fig. | we consider the same system as in Fig. || 
(one perfectly polarizing F layer, one partially polarizing 
layer), but as a function of the conductance of the normal 
layer for angles 9 close to 0. We find that switching 
the two layers has a drastic effect on the torque, even 
at a qualitative level. Interestingly, in the case where 
F a is the nearly perfect layer (dashed line in Fig. |]), a 
maximum of the torque is found for ^ < 1, i.e., in that 
case, a dirty metal spacer would give a higher torque (per 
unit of current) than a clean one. 

At this stage, it is interesting to compare our theory 
to that of Ref || In this work the WKB approxima- 
tion was used, and the electrons at the FN interfaces are 
either totally transmitted or reflected. Fori non-perfect 
polarizers, only a fraction of the channelsE3 act as per- 
fect filters while the others perfectly transmit both the 
minority and majority spins. However, this situation is 
different from having non perfect transmission probabil- 



ities T 



I' 



Tj per channel. 



In particular, having (Tj)iV 
channels that do not filter and (1 — (T|))7V perfect filters 
is not equivalent to N channels that all partly transmit 
the minority spins with probability (Tj). This situation 
is illustrated in Fig |7|. T-bjcJatter scenario is supported 
by ab-initio calculations .EjEj Moreover, for a disordered 
normal-metal spacer, multiple scattering from impurities 
mixes all channels and the notion of two type of channels 
become superfluous. In that case, the torque is described 
by Eq. (p3|) in all cases. The torque found in the second 
case can be significantly smaller than under the assump- 
tion of Ref ||. 



FIG. 6. Derivative of the torque with respect to 9 at 6 — 
as a function of #n, in unit of for the case where Fj, is 
a nearly perfect polarizer (|itt| 2 ~ 0.999, \tbi\ 2 = 0.001) and 
F a is not (|t aT | 2 = 0.3 \t al \ 2 = 0.01) (solid line), and for the 
opposite setup, F a is a nearly perfect polarizer (|t f j 2 = 0.999, 
\t al \ 2 = 0.001) while F 6 is not (\t b -j\ 2 = 0.3 \t bl \ 2 = 0.01) 
(dashed line). 



We can also compare our model to the wort of Berger.H 
While the theories of Berger and Slonczewsfl have much 
in common, Berger does invoke inelastic spin-flip scatter- 
ing in a way that Slonczewski does not. (Slonczewski's 
theory utilizes only spin-filtering, without spin-flip scat- 
tering.) This effect can in principle be treated in our 
model, by including the off-diagonal spin-flip reflection 
and transmission amplitudes that we have thus far ne- 
glected. We shall comment on some of the implications 
in the conclusion. We suspect that the differing treat- 
ments of this aspect of the physics may explain why Slon- 
czewski and Berger predict slightly different forms for the 
current-induced torques. 

In our theory, the scattering matrices of the ferromag- 
netic layers appear as free input parameters. However, it 
is in principle possible to calculate them from first princi- 
ple calculations for specific materials. Such an approach 
has been taken in Ref. and the results can be used 

to give some estimates of torques that can be expected 
in realistic systems. In Fig. |^, we compare the Co-Cu-Co 
system considered in the experiment of Ref. fj] with the 
Fe-Cr-Fe system. In the latter, the minority spins have 
a larger transmission probability than the majority ones, 
explaining the opposite sign of the torque. 



VI. CONCLUSION 

We have developed a theory for the spin-transfer- 
induced torques on the magnetic moments of a 
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FIG. 7. Torque per unit current as a function of 6. The 
solid line shows the case where the minority spin are trans- 
mitted with probability (Tj_) = 20% for all the channels. The 
dashed line shows the case where the minority spins are trans- 
mitted with probability (Tj) = 1 for 20% of the channels, 
otherwise (theory of Ref. H). The majority spin are totally 
transmitted in both cases. Inset: same system, dr^/IdO at 
6 = as a function of {T±} for the two different models. 
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FIG. 8. Torque per unit current as a function of 9 for 
two different realistic systems. The solid line shows the 
Co-Cu-Co trilayer (|t aT | 2 = 0.73, |i a i| 2 = 0.49, |t 6T | 2 = 0.68, 
\t H \ 2 = 0.29). The dashed line shows the Fe-Cr-Fe trilayer 
(|t aT | 2 = 0.48, \t al \ 2 = 0.59, |t 6T | 2 = 0.30, \t bL \ 2 = 0.50). In 
both cases = N has been assumed. The thick layer is 
assumed to be semi-infinite while for the thin layer, only the 
interface properties have been taken into account. Numerical 
values are obtained from Ref. El[ 



ferromagnet-normal-ferromagnet FNF trilayer system 
caused by a flowing current. Our theory deals with the 
effects of multiple scattering between the layers using the 
scattering matrices of the ferromagnet-normal metal in- 
terfaces as input parameters. We consider both the cases 
of a diffusive and ballistic normal metal spacer. Remark- 
ably, in the diffusive case, the high-dimensional scatter- 
ing matrices of the FN interfaces only appear through the 
reduced 4x4 tensor products of Eq.(J19|) which greatly 
reduces the number of degrees of freedom of the the- 
ory (see also Ref. [r|). This reduction of the number of 
degrees of freedom allows us to make qualitative predic- 
tions about the role of the interface transparency, normal 
metal resistance etc., without detailed knowledge of the 
microscopic details of the system. However, for quantita- 
tive predictions, inclusion of the microscopic parameters 
in our theory, e.g. from ab initio calculationsllinij is still 
needed. 

Having a complete theoretical description of the 
current-induced switching of magnetic domains in FN 
multilayers as a final goal, the theory here can be re- 
garded as being an intermediate step. On the one hand, 
microscopic input is needed for the scattering matri- 
ces of the FN interfaces, as explained above. On the 
other hand, the output of our theory, the current-induced 
torques, needs to be combined with restoring (hysteretic) 
forces in a more phenomenological theory that describes 
the dynamics of the magnetic moments. Such a theory in- 
volves anisotropy forces and information about the mech- 



anism by which the torque is exerted (spin wave excita- 
tion, locaLfijrhange field) - issues which are still subject 
of debate.Erd'Oo. 

In this paper, we have focused on the effects of "spin fil- 
tering" as the mechanism for current-induced torque, i.e., 
the difference in the transmission and reflection probabil- 
ities for electrons with spins parallel and antiparallel to 
the moments of the ferromagnetic layers (the diagonal 
terms in the matrices for the reflection and transmission 
amplitudes, Eq.([|).) A different source of spin-dependent 
scattering, which we have not considered in detail, but 
which is included in our formalism, is that of spin-flip 
scattering - the off-diagonal terms in Eq.(|). Its effect 
can be twofold. In the normal spacer, it would decrease 
the effective polarization, and therefore the torque. How- 
ever, in the ferromagnet, the rate of spin-flip scattering 
might be asymmetric with respect to minority and ma- 
jority spins, and therefore spin-flip scattering may also be 
an additional source of torque. As the number of degrees 
of freedom involved is much larger than for spin filtering 
only, a realistic model for the scattering matrices in the 
ferromagnets would be a necessary starting point for a 
theory that would include the effect of spin-flip scatter- 
ing. We leave such a theory for future work. 

ACKNOWLEDGMENTS 

We thank A. Brataas and G. Bauer for drawing our 
attention to Ref. O. 



9 



APPENDIX 

A. Derivation of Eq.(§|) 

In this appendix, we describe the calculation of Eq. ( pjj ) 
step by step. 

First, we substitute the expression @ for Q into (jl4|), 
and then formally expand the resulting equation in pow- 
ers of the reflection matrices r ai rj,, r n and r' n . Using the 
polar decomposition Eq. (|l7]) for the reflection and trans- 
mission matrices r„, i„, t' n and r' n of the normal layer, we 
get a sum of many terms, each of which is of a form where 
contributions from N are alternated with those of F a and 
Fb. Writing spin indices explicitly (summation over re- 
peated indices is implied), we can write those terms as, 



Tr 



N t l0 (A jk aB kl (3...r,C sm ) (Fl m J ...5^ D 



(36) 



where A, B, C, D, E, F e 
to the layer F Q and F^ while 
{uiVTv,uVl-Tu', 



{r a ,t a ,r b ,r' b ...} refer 
a, (3, 7, S, 77, uj € 
} refer to the normal layer. 



We are now ready to do the average of Eq. ((3^) over the 
matrices u, u' ', v and v' using the diagrammatic technique 
of Ref. (In leading order in TV these integrals reduce 
to the application of wick theorem.) Doing so, each of 
the a, (3, ... has to be put in correspondence with one of 
the y',8>, etc. To leading order in N, only the ladder 
diagram survives, in which a = 7, (3 = 6, rj = u),... and 
hence, A = D, B = E, C = F,... Thus, after averaging, 
we get terms like: 



Tr 



JY 



N ' 



-Tr 



N 



(37) 



where a, b, c, ... stands for either -^-TrT or j^Tr(l — T). 

To leading order in N, the average over T can now be 
done by simply replacing a, b, c, ... by their average value 
g-^/N or 1 — g-^/N where <?n is the average conductance 
(per spin) of the normal layer, in unit of e /h. 

Finally, denoting A = and /_t = (k,p), let us now 
introduce 4x4 matrices A, B, C,... that are defined as: 



AikA 



and S is defined as 



(38) 



(39) 



In term of these new matrices, eq.(37) now reads as a 
simple matrix product: 



T14ZA&BP.. .fjC, 



with a, $, fj,...e {gn/N,l - g N /N}. Equation @ is 
formally equal to the expansion of f2 (see Eq.(^6|)) except 
that we are now dealing with "hat" matrices. Therefore, 
we can now resum all the terms of the expansion and get 
Eq. 



B. Ballistic normal layer: a pedestrian approach 

If N is very clean, and the interfaces are very flat, it is 
reasonable to assume that the electrons propagate ballis- 
tically inside the normal layer. The different modes will 
not be mixed in that case, and the electron wavefunc- 
tion only picks up a phase factor e lkiL where L is the 
width of N and ki the momentum of channel i. For a 
sufficiently thick normal layer (i.e. L ^> Xp), small fluc- 
tuations of k{ lead to an arbitrary change in the phase 
factor, and it is justified to consider e lkiL as a random 
phase and to average over it. This is different from the 
case of a disordered metal spacer, where the average in- 
volves unitary matrices u, u' ... that mix the channels, cf 



Eq.(T' 



In the case where r. 



are proportional 



to the identity matrix (i.e. the reflection amplitudes do 
not depend on the channel), the ballistic model reduces 
to the disordered model of Eq. (^J) f° r 9n = N. 

The reflection matrices of N being zero, the matrix 
reads: 



JkiL 



1 



1 



(41) 



Neglecting spin-flip scattering, denoting z — e 
choosing r all = r ah r a22 = r al ,... where r aT , r al 
diagonal matrices, one gets after some algebra: 



2ikiL anc | 

are 



dV ( ' 



ev A(z) 
Tr Re — — 

4tt z\D(z)\ 2 



(42) 



where A(z) and D(z) stand for: 

A(z)=(l-t' b / b \-r' M r' b \) 

( |f oT | 2 (l-^r a i)(^-^) 
- \t' ai \\l-zr H r^){z-r'^)) 



(43) 



D(z) 



cos2 o(M r M + r ai r 'b\) 



+ sin 2 -(r oi r^ T + r al r' H ) 
+ z 2 r ai r a] r' bl r' bV 



(44) 



A similar formula can be written for the conductance 
9(z): 



= £! T b{z) 

9 h 1 z\D(z)\ 2 '' 



(45) 



(40) with: 
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B{z) = \t' a] 
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2 cos 2 


2 (1 ~ 


zr air' H )(z 


~ r al r bl> 


+ K T 


% 


2 sin 2 


I' 1 - 

2 




* \ 

— r oi r 6t/ 


+ Ki 


% 


2 sin 2 


fo- 


zr a ^r' bl ){z 


* \ 


+ Ki 


% 


2 cos 2 


ld- 


zr a]r' b] ){z 


~~ r a^ r bl) 



Taking the average over the phases now amounts to cou- 
tour integration for z: 



(47) 



where the integration is done along the unit circle. The 
result is then given by the sum of the poles that are inside 
the unit circle. The two poles of D(z) are outside the unit 
circle, while the two poles z\ and Z2 of z 2 D (—) are inside 
the circle. They are given by: 

Zi = Tjcos 2 ^(r aT r£ T + r al r' H ) + sin 2 ~(r aJ r£ t + r a] r' b[ ) 



cos4 o( r «T r 6T - r ai r 'b\) 



2 cos 2 - sin 2 -(r' bl r' b1 (r a1 



+ sin 4 -(r al rL - r a1 r' b ,) 2 



fair ^ (r; 



^)) 

(48) 



The averaged torque and conductance are then simply 
given by 



and 



Wo' 
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4-7T 21 — Zi 



-Tr 



h Zl — Z2 



-Tr 



A(zi) A(z 2 ) 

I>(2l) ^(Z 2 ) 

gfa) B(za) 



(49) 



(50) 



In the case where all the channels are not identical, these 
results can be generalized by introducing a k dependence 
of the different transmission/reflection amplitudes. 
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